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Abstract 



We report the results of a study of color coherence effects in pp collisions 
based on data collected by the D0 detector during the 1994-1995 run of 
the Fermilab Tevatron Collider, at a center of mass energy ^/s = 1.8 TeV. 
Initial-to-final state color interference effects are studied by examining particle 
distribution patterns in events with a W boson and at least one jet. The 
data are compared to Monte Carlo simulations with different color coherence 
implementations and to an analytic modified-leading-logarithm perturbative 
calculation based on the local parton-hadron duality hypothesis. 
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Color coherence phenomena in the final state have been studied since the early 1980's 
in e + e~ annihilations [HH^] and are very well established. The study of coherence effects in 
hadron-hadron collisions is considerably more subtle than those in e + e~ annihilations due 
to the presence of colored constituents in both the initial and final states. In this paper 
we report the first results on initial-to-final state color interference effects in pp interactions 
using jy+jets events. 

Coherence phenomena are an intrinsic property of any gauge theory. In quantum chro- 
modynamics (QCD), color coherence phenomena can be instructively separated into two 
regions: intrajet and interjet coherence |7],[J. Intrajet coherence deals with coherent effects 
in partonic cascades, resulting on average in the angular ordering (AO) of the sequential 
parton branches, which give rise to the depletion of soft particle production (the so called 
"hump-backed plateau") inside jets pHIIf . Interjet coherence is responsible for the string 
|T^] or drag ]13| effect first observed in the final state products of e + e~ annihilations. It 
deals with the angular structure of soft particle flows when three or more energetic partons 
are involved in the hard process. In this case, the overall structure of particle angular dis- 
tributions is governed by the underlying color dynamics of the hard scattering processes at 
short distances. 

Perturbative quantum chromodynamics (pQCD) calculations have been used to describe 
the production of jet final states. However, descriptions of the characteristic particle struc- 
ture of high energy hard collisions still rely on phenomenological models to explain how the 
partonic cascade evolves into final state hadrons. Within Monte Carlo (MC) simulations 
incorporating such models, the primary partons from the hard scatter evolve into jets of par- 
tons via gluon and quark emission according to pQCD. This process continues until a cut-off 
transverse-momentum scale (Qq ~ 1 GeV/c) is reached. After this phase, non-perturbative 
processes take over, which "cluster" the final partons into color singlet hadronic states via 
a mechanism described by phenomenological fragmentation models, like the Lund "string" 
rj| or the "cluster" |T5| fragmentation models. These simulations usually involve many a 



priori unknown parameters that need to be tuned to the data. 

A different and purely analytical approach giving quantitative predictions of hadronic 
spectra is based on the concept of Local Parton Hadron Duality (LPHD) ||. The key 
assumption of this hypothesis is that the particle yield is described by a parton cascade in 
which the conversion of partons into hadrons occurs at a low virtuality scale, of the order 
of hadronic masses (Qq ~ 200 MeV/c 2 ) and independent of the scale of the primary hard 
process, and involves only low momentum transfers. It is assumed that the results obtained 
for partons apply to hadrons as well in an inclusive and average sense. Within the LPHD 
approach, resummed pQCD calculations for the parton cascade have been carried out in the 



simplest case (high energy limit) in the Double Logarithmic Approximation (DLA) |]T0| , p~6 



and in the Modified Leading Logarithmic Approximation (MLLA) P, |l7l , |T8f , which includes 
higher order terms of relative order ^foT s (e.g., finite energy corrections). These higher order 
terms are essential for quantitative agreement with data at present energies |7],[n| . 



The AO approximation is an important consequence of color coherence. It results in the 
suppression of soft gluon radiation in partonic cascades in certain regions of phase space. For 
the case of outgoing partons, AO requires that the emission angles of soft gluons decrease 
monotonically as the partonic cascade evolves away from the hard process. MC simulations 
including coherence effects probabilistically by means of AO are available for both initial 
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and final state evolutions. (Parton shower event generators incorporate AO effects in the 
initial state as the time reversed process of the outgoing partonic cascade, i.e. the emission 
angles increase for the incoming partons as the process develops from the initial hadrons to 
the hard subprocess.) AO is an element of the DLA and MLLA analytic pQCD calculations, 
which provides the probabilistic interpretation of soft-gluon cascades. In fact, beyond the 
MLLA a probabilistic picture of the parton cascade evolution is not feasible due to 1/N%- 
suppressed (where N c is the number of colors) soft interference contributions that appear in 
the higher-order calculations ]8|,|T7i|. 

Both the CDF [2U| and D0 [2~I] Collaborations have measured spatial correlations be- 
tween the softer third jet and the second leading-^ jet in pp — > 3 jets + X events to 
explore the initial-to-final state coherence effects in pp interactions. The extraction of the 
color coherence signal in these measurements relies on comparisons of data distributions 
to MC simulations with and without coherence effects. In the analysis described here, the 
coherence signal in the data is extracted in a more direct way by comparing the soft particle 
angular distributions around the colorless W boson and opposing jet in the same event. 



The D0 detector is described in detail elsewhere [52]. This analysis uses the tracking 
system and the uranium/liquid- argon sampling calorimeter. The D0 calorimeter has a 
transverse granularity of Arj x A0 = 0.1 x 0.1 forming projective towers, where r\ is the 
pseudorapidity (77 = — ln[tan(6 l /2)], 9 is the polar angle with respect to the proton beam), 
and is the azimuthal angle. It has hermetic coverage for \rj\ < 4.1 with fractional transverse 



energy Et resolution of ~ 80%/ J -Sr(GeV) for jets and fractional energy resolution of 



15%/^£(GeV) for electrons. 

The data sample for this analysis P31, representing an integrated luminosity of 85 pb _1 , 



was collected during the 1994-1995 Tevatron Collider run. Events from W — > e + v 
decays were collected with a trigger that required a minimum missing transverse en- 
ergy ($ T ) of 15 GeV and an isolated electromagnetic (EM) cluster with transverse en- 
ergy Et > 20 GeV. The offline kinematic requirements imposed on this sample were 
$r > 25 GeV, E T > 25 GeV, \rj e \ < 1.1, and 40 GeV/c 2 < M T (e,$ T ) < n GeV/c 2 , 



where MT(e,$ T ) = y 2pr(e)pr(i / )(l ~ cosA0)/c 2 is the W boson transverse mass and A0 
is the azimuthal separation between the electron and neutrino. The transverse momentum 
of the neutrino, Pt{v), was calculated using the calorimetric measurement of the ~$ T in the 
event. 

Events were required to have one electron cluster passing four quality criteria based on 
shower profile and tracking information: (i) the ratio of the EM energy to the total shower 
energy had to be greater than 0.95, (ii) the position of the calorimeter energy deposition of 
the electron had to match with a track found in the drift chambers, (iii) the lateral and lon- 
gitudinal shape of the energy cluster had to be consistent with those of an electron, and (iv) 
the electron had to be isolated from other energy deposits in the calorimeter with isolation 
fraction / iso < 0.1. The isolation fraction is defined as /i SO = [£7(0.4) — £ , em(0.2)]/£'em(0.2), 
where E(R cone ) (-EEM(-Rcone)) is the total (electromagnetic) energy within a cone of radius 



-Rcone = y (A77) 2 + (A0) 2 centered around the electron. 

Jets in the events were reconstructed offline using an iterative fixed-cone clustering algo- 
rithm with cone radius -R con e = 0.7. Spurious jets from isolated noisy calorimeter cells and 
accelerator losses were eliminated by loose cuts on the jet shape. 
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Events were required to have a measured vertex with longitudinal position within 20 cm 
of the detector center to preserve the projective geometry. Since multiple interactions (more 
than one proton-antiproton interaction in the same bunch crossing) are expected to increase 
the global energy level in the calorimeter affecting the color coherence signal, we retained 
only events with a single reconstructed vertex and additionally required the beam-beam 
hodoscope timing information to be consistent with a single interaction. Finally, events 
were eliminated when there was significant pileup energy in the calorimeter around the 
region where the Tevatron Main Ring passes through the D0 detector. 

We study color coherence in H^+jets events by comparing the distributions of soft par- 
ticles around the W boson and opposing leading-.E r jet. Since the W boson is a colorless 
object, it should not contribute to the production of secondary particles, thereby providing 
a template against which the pattern around the jet may be compared. This comparison 
reduces the sensitivity to global detector and underlying event biases that may be present 
in the vicinity of both the W boson and the jet. 

The W boson was reconstructed from its electron and neutrino decay products, resulting 
in a twofold ambiguity in the W boson rapidity yw (due to the corresponding ambiguity in 
the neutrino longitudinal momentum pz). MC studies have shown that the smaller \yw\ is 
closer to the true W boson rapidity approximately 2/3 of the time, so this is the solution 
chosen. This choice was also made in the MC W boson reconstruction to retain consistency 
in the comparison with the data. 

Once the W boson direction has been determined, the opposing jet was identified by 
selecting the leading-^ jet in the azimuthal hemisphere opposite to the W boson. Annular 
regions are drawn around both the W boson and the tagged jet in (77, 0) space. The angular 
distributions of calorimeter towers with E?> 250 MeV are measured in these annular regions 
using the polar variables R = ^(Ar/) 2 + (A0) 2 and (3 X = tan^^g^); where X = W 
or jet, Ar] W = r/tower - yw, &Viet = ^tower - Viet, and A<fix = 0towcr - (fix, in a search disk of 
0.7 < R < 1.5 (Fig. [l]). For events in which the jet and the W boson annuli overlapped, all 
calorimeter towers in the shared region were assigned to the nearest object. 

We define /3w(jet) = to point along the beam direction nearest to the W bo- 
son (jet). Calorimeter towers which correspond to partially instrumented calorimeter re- 
gions (regions between the central and end-cap calorimeter cryostats with pseudorapidities 
1.1 < |r/ towcr | < 1.4) were not included in the (3 distributions for either the data or MC sim- 
ulations. Calorimeter cells that belong to the electron cluster in a cone of R — 0.3 from the 
centroid center were eliminated. We study the interference effects in regions \Vjet\ < 0.7 and 
\yw\ < 0.7, requiring the tagged jet to have Et > 10 GeV and the W boson px > 10 GeV/c. 
After application of all selection criteria 390 events remain. 

The measured angular distributions are compared to the predictions of PYTHIA 5.7 |24j 
parton shower event generator with different levels of color coherence effects and to an ana- 
lytic pQCD calculation of Khoze and Stirling E3| based on MLLA and LPHD. The PYTHIA 



MC sample was processed through a full GEANT-based detector simulation f2q| . To best 
model the calorimeter noise effects in the MC simulations, we overlaid noise contributions 
for each calorimeter cell from the data. The generated events were subsequently processed 
using the same criteria employed for analyzing the data. 

PYTHIA incorporates initial and final state color interference effects by means of the AO 
approximation of the parton cascades. After the perturbative phase it employs the Lund 
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string fragmentation (SF) model (or an independent fragmentation IF model) as the phe- 
nomenological model to describe the non-perturbative hadronization process. The SF model 
has been supported by the observations of color coherence phenomena in e + e~ annihilations. 
In PYTHIA, the AO constraint can be turned off. When both AO and SF are implemented, 
PYTHIA accounts for color coherence effects at both the perturbative and non-perturbative 
levels. Turning off AO removes the perturbative contribution, and using IF eliminates the 
non-perturbative component. 

We check whether the MC simulations describe the event characteristics in the data 
using distributions of electron Et, the event $ T , and the azimuthal and rapidity separation 
of the W boson and the tagged jet. These are compared to the PYTHIA simulations with 
full coherence effects in Fig. |2|. For all distributions, PYTHIA is in good agreement with the 
data. 

Figure ^|(a) shows the measured angular distributions of the number of towers above 
threshold around the jet and around the W boson. The Et of each calorimeter cell (each 
calorimeter tower is constructed from many calorimeter cells following the projective geom- 
etry of the detector with respect to the detector center) was corrected for offsets due to 
noise, zero suppression, and energy pileup effects ||23|| . A prominent feature of both curves is 
a strong peaking around n/2 < (3 < 37r/4. This is because the shapes of the (3 distributions 
are sensitive both to process dynamics and to phase space effects resulting from our event 
selection criteria and the calorimeter tower Et threshold. 

Figure 0(b) shows the ratio of the tower multiplicity around the jet to the tower multiplic- 
ity around the W boson as a function of (3. The errors include only statistical uncertainties, 
which are the dominant source of uncertainty. The data show that the tower multiplicity 
around the jet is greater than that around the W boson and the excess is enhanced in the 
event plane (i.e., the plane defined by the directions of the W boson or jet and the beam 
axis: (3 = 0, 7r) and minimized in the transverse plane (j3 = 7r/2). This is consistent with the 
expectation from initial-to-final state color interference that there is an enhancement of soft 
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FIG. 2. Comparisons of data (points) with PYTHIA events with AO and SF (histogram) for 
(a) electron Et, (b) event $ T , (c) azimuthal separation of the W boson and tagged jet, and (d) 
rapidity separation of the W boson and tagged jet. 



particle production around the tagged jet in the event plane relative to the transverse plane 
when compared with the particle production around the W boson. It is also in agreement 



with our published multijet analysis results [21 



The ratio of the tower multiplicity around the jet to the tower multiplicity around the 
W boson is shown in Fig |] for the data and the MC predictions as a function of (3. All 
predictions have been normalized to the integral of the data /3-distribution. PYTHIA with 
AO and SF is in good agreement with the iy+jets data. PYTHIA with AO off and SF 
agrees less well, and PYTHIA with AO off and IF does not reproduce the data. The shape of 
the analytic prediction based on MLLA and LPHD is consistent with the data as shown in 
Fig. |](d). We note that the analytic calculation was performed assuming that the W boson 
and the outgoing parton were produced centrally and does not include any underlying event 
or detector simulation effects. 

To measure the color coherence signal, we construct the variable R sig as the jet/W boson 
tower multiplicity ratio of the event plane {(3 = 0,ir) to the transverse plane {(3 = tt/2). 



We define R sig = f»|, where R{(3 = 0,tt) = gl^Z , R((3 = tt/2) = |£, 

and N^^2) ls the number of towers above threshold for the i th bin of the W boson (jet) 
P distribution. i? sig is expected to be near unity in the absence of color coherence effects. 
In addition R sig is insensitive to the overall normalization of the individual distributions, 
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FIG. 3. (a) Calorimeter tower multiplicity around the jet and W boson as a function of 0. (b) 
Ratio of the tower multiplicity around the jet to the tower multiplicity around the W boson as a 
function of (3. 



and MC studies have shown that it is relatively insensitive to detector smearing effects. 
Figure |5] compares the _R S ig variable for the data to the various PYTHIA predictions and to 
the MLLA+LPHD calculation. Clearly the value of R s \ g for the data deviates from unity 
in agreement with PYTHIA with AO on and SF, and in disagreement with AO off and 
SF or AO off and IF. These comparisons imply that for the process under study, string 
fragmentation alone cannot describe the effects seen in the data. The AO approximation is 
an element of parton-shower event generators that needs to be included if color coherence 
effects are to be modeled successfully. Finally, the analytic prediction by Khoze and Stirling 
is consistent with the data, thus giving additional evidence supporting the validity of the 
LPHD hypothesis. 

The dominant source of uncertainty on the data f3 distributions is statistical due to the 
limited event sample. The statistical error for R sig is 7%. Since we report ratios of event 
distributions any possible uncertainty on quantities that affect the overall rate of events 
is minimized. Sources of systematic uncertainty arise from background contamination to 
iy+jets events, uncertainties in the calorimeter channel-to-channel offset correction, multi- 
ple pp interactions, and uncertainties associated with the calorimeter tower threshold. 
The primary background to VT+jets events is dijet production in which one of the jets mim- 
ics the characteristics of an electron. The $ T in such events typically arises from shower 
fluctuations or calorimeter imperfections. For our selection criteria, the estimated back- 
ground level is about 5% resulting in a 1-2% uncertainty on the coherence signal. The effect 
on the -R S ig variable due to uncertainties in the offset correction were found to be at the 
1% level. To evaluate possible effects on the signal coming from residual multiple interac- 
tion contamination, we examined the dependence of i? S i g as a function of luminosity. No 
systematic dependence of the signal variable was observed for our event sample. 

The dependence of the signal variable i? S i g on the calorimeter tower threshold was 
studied by varying the threshold from 200 to 350 MeV. Figure ^| shows how the R sig variable 
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FIG. 4. Comparison of the Jet/VF boson tower multiplicity ratio from data to (a) PYTHIA with 
AO on and SF, (b) AO off and SF, (c) AO off and IF, and to the (d) MLLA+LPHD predictions. 
The predictions have been normalized to the data. All uncertainties are statistical only. 



varies as a function of the tower threshold for the data and for the three PYTHIA predictions. 
Although there seems to be a dependence of i? S i g with the tower Et threshold for both the 
data and the simulations, for all tower thresholds examined the PYTHIA predictions with 
AO on and with SF resemble the data best. 

In summary, color coherence effects in pp interactions have been observed and studied 
by the D0 Collaboration. We have presented the first results on color coherence effects 
in H/+jets events. The data show an enhancement of soft particle radiation around the 
jet in the event plane with respect to the transverse plane, consistent with color coherence 
as implemented in the PYTHIA parton shower event generator, which includes the angular 
ordering approximation and string fragmentation. In addition, the relative amount of en- 
hancement is consistent with an analytic perturbative QCD calculation based on modified 
leading logarithmic approximation and local parton-hadron duality. 
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Sjostrand for helping us with the color coherence implementations in PYTHIA. We thank the 
Fermilab and collaborating institution staffs for contributions to this work, and acknowledge 
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sariat a L'Energie Atomique (France), Ministry for Science and Technology and Ministry 
for Atomic Energy (Russia), CAPES and CNPq (Brazil), Departments of Atomic Energy 
and Science and Education (India), Colciencias (Colombia), CONACyT (Mexico), Ministry 
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